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Abstract 

The paper intends to propose an algorithm which could identify a general class of pdes 
describing dynamical systems with similar symmetries. The way that will be followed starts 
from a given group of symmetries, the determination of the invariants and, then, of the 
compatible equations of evolution. The algorithm will be exemplified by two classes of 
equations which describe the Fokker-Planck model and the "backward" Kolgomorov one. 

1 Introduction 

Lie theory of symmetry groups for differential equations is one of the most important methods 
for studying nonlinear problems which appear in physics or in other fields of applied mathe- 
matics. The main idea of Lie's theory is the investigation of integrability starting from the 
invariance of the equation under some linear transformations of independent and dependent 
variables, transformations which define the so called Lie group of symmetries. There are dif- 
ferent types of symmetries which can be identified for different differential equations: point-like 
symmetries, contact symmetries, classical, generalized or non-local symmetries. Many studies 
have been devoted to this topic, some of them being [HEIIH]. The existence of the Lie symmetry 
generators for differential equations often allows the reduction of those equations to simpler 
ones. The similarity reduction method for example may be a way to follow for transforming 
partial derivative equations (pdes) defined in 2D space in ordinary differential equations (odes). 

Usually, the direct symmetry problem of evolutionary equations is considered. It consists 
in determining the symmetries of a given evolutionary equation. The aim of this paper is 
to investigate the inverse problem: what is the largest class of evolutionary equations which 
are equivalent from the point of view of their symmetries. We propose an algorithm which 
allows finding a general class of pdes describing dynamical systems with similar symmetries. 
The problem will be effectively formulated in the section 2 of the paper for an equation of 
the form ut{x,t) = F{t,x,u,Ux,U2x)- A general algorithm for identifying the most general 
equation with a given symmetry group is presented. It is also pointed out how this pde defined 
in (1 + 1) dimensions can be reduced to an equivalent ode by using the similarity reduction 
procedure. The general results will be particularized in the section 3 of the paper for two 
important models represented by the Fokker-Planck equation and "backward" Kolgomorov one. 
Generalized Fokker-Planck equations can be derived from generalized linear non-equilibrium 



thermodynamics [3]. It was already shown that the Fokker-Planck equation in 2D could be 
transformed, by appropiate change of coordinates, in heat equation or Schrodinger equation [5]. 
As far as Kolmogorov equation is concerned, it has been considered even in infinite dimensions 
and interesting applications to stochastic generalized Burgers equations have been noted [6]. 
Some concluding remarks will end the paper. 

2 The general setting of the problem 

Let us consider a general dynamical system described in a (1 + 1) dimensional space {x, t) by a 
second order partial derivative equation of the form: 

Ut = A{x,t,u)u2x + B{x,t^u)Ux + C{x,t,u)u, A (1) 

The coefHcient functions A{x,t,u), B{x,t,u), C{x,t,u) are arbitrary functions. Any equa- 
tion Ai of the form ([T|) is invariant under the infinitesimal transformations: 

x = X + ^{t,x,u)£ + t = t + ip{t,x,u)£ + u = t + r]{t,x,u)e + ... 
if and only if it satisfies [7] : 

t/^")(Ai) 1a,=o=0 (2) 

where t/*^"^ is the n— th extension of the symmetry operator U. In the case of the equation ([1]) 
defined in (1 + 1) dimensions, the operator U, also known as Lie operator, has the general form: 

Odd 
U{x,t,u) = ip{x,t,u)— +^{x,t,u)— + r]{x,t,u)— (3) 

As in our case the equation ([1]) is a second order differential equation, the invariance condition 
([2]) win be in fact: 

C/*-^-' [ut — A{x, t, u)u2x — B{x, t, u)ux — C(x, t, u)u] = (4) 

where C/*-^^ is the second order extension of the Lie symmetry operator U . A concrete computa- 
tion leads to the equation: 



= {^At + iAx + 'nAu)u2x + {^Bt + iBx + 'nBu)ux + ^uCt + iCxU+ (5) 
+Cri + r]CuU + Bt]'' - -rf + Aif"^ 

By substituting in d^D the general expressions given in [7] for rj^ , rf, rf'^ and asking for the 
vanishing of the coefficients of each monomial in the derivatives of u{t,x), one obtains the 
following differential system: 

= 0; = 0; = 0; r/2„ = 0; -ipAt - ^Ax - r^A^ - Aip, + 2Ai^ = 0; 
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- ipBt - iB^ - vBu + B^, -Ct- Bif^ - 2A7],^ + 71^2^ = 
-ifCtu — £,CxU — Ct] — r]CuU — Br]^ + % + Cij^u — (ffCu — vlr/2x = 



(6) 



It is important to note that this system can be seen as a general system with the unknown 
functions A{x,t,u), B{x,t,u), C{x,t,u), ip{x,t,u), r]{x,t,u).Two completely different 

situations can be considered: 

- one fix the equation ([1]) by choosing concrete expressions for A{x, t, u), B{x, t, u), C{x, t, u) 
and one determines the form of the symmetries, that is ip{x,t,u), ^(x,t, n), r]{x,t,u); 

- one looks for the class of equations which have a given form of the symmetry, which 
means to solve (0) with given ip{x,t,u), ^{x,t,u), r]{x,t,u) for the unknown functions A{x,t,u), 
B{x, t, u), C{x, t, u). 

Usually the first approach is considered. What we are doing here is to perform the second 
approach. We will start from the already known symmetries of the Fokker-Planck equation [1], 
an equation which belongs to the general class of equations ([TJ and describes various evolutionary 
processes in quantum optics, solid state physics or statistical physics. Let us therefore consider 
for the Lie symmetry operator ([3]) the following form: 



„)=.-| + .e^.|-.V'.|; (7) 
With these choices for ip(x,t,u), ^{x,t,u), r]{x,t,u), the system ([6]) reduces to the equations: 



At + xAx — X AuU = 
-Bt - xB^ + x^uBu - B -2x + 4Ax = 
CfU + Cxux - CuU^x^ - 2Bxu + 2x^u + 2Cu - 2Au = 



(8) 



It is a system with the unknown functions A{x,t,u), B{x,t,u), C{x,t,u) which leads to the 
general solution of the form: 



A = f (^t — Inx, ue 



x2/2 



B 
C 



-1 + 2/ \^t-lnx,ue 
—l + f(^t — \nx, ue 
h (t — Inx, ue^'^l'^ 



X + 



g it — \nx, ue 



(9) 



+ / ~ ln3;,ne' 



g\t — \nx,ue 
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The solution ([9|) is given in terms of 3 arbitrary functions f{u,x,t), g{u,x,t), h{u,x,t). In fact, 
two types of combinations of the variables {u,x,t) have to be considered: 



z= (t - In x) , = ne"^^/^ 



(10) 
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The significance of these variables becomes clear if one considers the characteristic equations 
associated with the symmetry operator ([7]): 

dt dx du ^^^^ 

By integrating the equations ()lip the following similarity variables are obtained: 



It is simple now to notice the direct connection between (jlOp and (jl2p . that is the general solution 
^ is given in terms of similarity variables. Moreover, it is also important to note that, by using 
the similarity transformations ()12p . the class of the evolutionary systems which are described 
by the (1 + 1) equation ([T|) can be reduced to the \D systems with the evolution given by the 
following differential equation: 

z^'^^f [0(z),-lnz] + z^<7 [<A(z),-lnz]+</>(z)/i[0(z),-lnz]=O, ^ f,g,h (13) 

As a matter of fact, let us note that for / = 1, (7 = /i = 0, the equation (fT3|) becomes: 

d'^Mz) , , 

One obtains a wave type equation in one dimension with the solution of the form: 

0(z) = ClZ + C2 (15) 

where ci,C2 are arbitrary constants. 

As a conclusion of this section, let us mention that the solution ^ gives the most general 
class of 2D equations of the form ([T]) which admits a symmetry operator of the form ([T]) and, by 
that, accepts a similarity reduction to the ID equation ()13p . In the next section we will identify 
two concrete examples of systems which belong to this class of evolutionary equations. 



3 Applications 

3.1 The Fokker-Planck model 

Let us come back to the choice / = l,(7 = /i = 0in the equation (jlSp . choice already considered 
before. We have to note that for this case the solution ([9]) of the system ([8|) transforms into very 
simple expressions: 

A = l, B = x, C = 1 (16) 
The equation ([T]) takes, it too, a very simple form known as the (1 + 1) Fokker-Planck equation 

Ut = U2x + XUx + u (17) 
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One recovers so that the equation ()17p is one of the particular form of the equation ([T]) which 
accepts the symmetry ([7]). As we mentioned, this is an already known result and, in fact, 
it justifies the choice of the symmetry operator in the form ([7]). The equation ()17p is used in 
various fields of natural sciences, as physical chemistry or theoretical biology [9j. In the similarity 
variables ()12p it can be reduced to form ()14p and can be solved. Coming back to the original 
variables x,t and u = u{x,t) by inverting (fT2]l . the solution (fTSj) of the equation becomes: 

n(x, t) = cixe-*e-^-'/2 + cae'^''/^ (18) 



So, it was possible to obtain a solution of the equation (|T8|) as a simple case of the general 
results obtained for the equation ([1]). We have to mention that, by computational methods, one 
can obtain a general solution of the equation (jlSp in a very complicated form. Such results are 
mentioned in [10], |11] . 

3.2 The "backward" Kolmogorov model 

Let us consider now another particular case of solution in ([9]): f = 1, g = —l,h = Q. In this case 
the equation ([T]) takes the form of an evolutionary equation of the type "backward" Kolmogorov 

ut = U2x + {x (19) 

X 

As the whole class of equations belonging to ([1]) and as the Fokker-Planck equation, the equation 
(|19p accepts the Lie symmetry ([7]) and can be reduced to a one-dimensional ordinary differential 
equation of the form: 

,fp^_d_m = o (20) 
az^ dz 

The equation ()20p has a solution of the form: 

0(z) = C3 + C4Z^ (21) 

Using again ()12p one transforms ()2ip into a solution of ()19p of the form: 

w(x, t) = e-^'/2 ^ c4x2e-2*) (22) 



4 Concluding remarks 

The results of this paper can be synthesized as follow: (i) the class of (1 + 1) generalized 
evolutionary equations of type ([1]) can allow the same concrete Lie symmetry operator ([7]) if 
and only if the coefficient functions A{x,t,u), B{x,t,u),C{x,t,u) have the expressions Q; (ii) 
the knowledge of symmetry operators enabled to introduce, by integrating the characteristic 
equations (|lip . the so called similarity variables. Through this similarity approach, the general 
form for reduced odes (fT3]l associated to the analyzed equations ([1]) has been obtained; {in) 
the results are particularized for two dynamical equations, the Fokker-Planck model and the 
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"backward" Kolgomorov one. The general algorithm allows deriving simple solutions of the 
evolutionary equations which describe the two models. 
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